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We introduce a concept of asymptotic principal values which enables us to handle rigorously 
singular integrals of higher-order poles encountered in the computation of various quantities based 
on correlation functions of a vacuum. 

Several theorems on asymptotic principal values are proved and they are expected to become 
bases for investigating and developing some class of regularization methods for singular integrals. 

We make use of these theorems for analyzing mutual relations between some regularization meth- 
ods, including a method naturally derived from asymptotic principal values. It turns out that the 
concept of asymptotic principal values and the theorems for them are quite useful in this type of 
analysis, providing a suitable language to describe what is discarded and what is retained in each 
regularization method. 
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I. INTRODUCTION 

Physics of quantum vacuum fluctuations is one of 
the intriguing research topics expected to be developed 
through the interplay between theories, experiments and 
practical applications. 

Investigations of quantum vacuum fluctuations even 
stimulate the border area between physics and mathe- 
matics. As a typical example of this sort, we often en- 
counter singular integrals in computing several quantities 
based on correlation functions of a vacuum in question. 
The occurrence of singularity or divergence is often a sig- 
nal of surpassing the border of validity of a model by too 
much extrapolation. Furthermore, it could originate from 
deeper physical processes for which satisfactory consis- 
tent mathematics is still unavailable. How to handle sin- 
gular integrals can be then a challenging topic, requiring 
both mathematical analysis and physical considerations. 

Faced with singular integrals, we need to resort to some 
regularization method to get a finite result. The aim of 
this paper is to give an organized mathematical basis un- 
derlying some typical regularization methods and make 
clear their mutual relations. We introduce below a con- 
cept of asymptotic principal values which can be a key 
tool to analyze some class of regularization methods. We 
then prove several theorems on the asymptotic principal 
values useful for studying regularization procedures. 

There are still various uncertainties to clear up in reg- 
ularization methods, reflecting our lack of mathematical 
basis for handling infinities. In this situation, we cannot 
expect any universal regularization method, but we need 
to customize the method by try and error depending on 
the problem in question. It is far from the aim of this 
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paper to judge which method is better than the others. 
It just tries to present a concrete mathematical basis for 
further considerations and developments of better regu- 
larization methods. 

The organization of the paper is as follows. In Sec. II, 
we present one simple example where singular integrals 
of a higher-order pole emerge. The origin of singularities 
in this example is physically clear and we can get some 
idea on how these integrals should be regularized. In Sec. 
Ill, we introduce asymptotic principal values which de- 
scribe precisely how singular integrals behave, which can 
be useful for investigating and developing several regular- 
ization methods. We then prove several useful theorems 
on the asymptotic principal values. In Sec. IV, we an- 
alyze some typical regularization methods by means of 
theorems prepared in the previous section. Sec. V is 
devoted for a summary and several discussions. 

II. TYPICAL EXAMPLE OF SINGULAR 
INTEGRALS 

Let us consider as an example the measurement of 
the electromagnetic vacuum fluctuations in a half-space 
bounded by a perfectly reflecting infinite mirror. Re- 
cently the switching effect [1] and the smearing effect 
due to the quantum spread of a probe particle [2] have 
been analyzed by studying the measurement process of 
the Brownian particle released in this environment. 

We thus take a model introduced in Ref . [3] and reana- 
lyzed in Ref.'s [1] and [2]: Suppose that a flat, infinitely 
spreading mirror of perfect reflectivity is placed on the 
xy-plane (z = 0). Then let us investigate the quantum 
vacuum fluctuations of the electromagnetic field inside 
the half-space z > by releasing a classical charged 
probe-particle with mass m and charge e in the envi- 
ronment. We can estimate the quantum fluctuations of 
the vacuum through the velocity dispersions of the probe- 
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particle released in the environment. 

When the velocity of the probe-particle is much smaller 
than the light velocity c, the motion for the particle is 
described by 

dv j=> 

m— = eE(x,t) , (1) 

where E(x,t) is the electric field. 

Within the time-period when the particle does not 
move so much, Eq. (1) along with the initial condition 
v(0) = vo is solved approximately as 

e r* - 

v(t)~v + — / E{x,t')dt' . (2) 
m Jo 

For simplicity, let us consider only the "sudden- 
switching" case; the measurement is switched on 
abruptly, stably continued for r [sec] before switched off 
abruptly. It is mathematically described by a step-like 
switching function without any switching tails. The ve- 
locity dispersions of the particle, {Av 2 ) (i = x, y, z), are 
then given by 

{Av 2 (x,t)) = ii f dt' [ dt"(E l (x 1 t')E l {x,t")) R , 
m Jo Jo 

(3) 

by noting that (Ei{x, t)) R = 0: Here (Ei(x, t')Ei(x, t")) R 
(i = x, y, z) are the renormalized two-point correlation 
functions of the electric field (the suffix "R" is for "renor- 
malized"). Now {Ei(x,t')Ei(x,t")) R (i = x,y,z) are 
computed [4] as 

(E z (x,t')E z (S,t"))n = ^ (r2 _ 1 (2z) 2 ) 2 W 
(E x (x,t')E x (x,t")) R = (E y (x,t')E v (x,t"))R 

7T 2 (T 2 - (2z) 2 ) 3 ' 1 ; 

where T ;= t' — t" . (We set c = h = 1 hereafter through- 
out the paper.) 

It is obvious that the integral in Eq.(3) is regular when 
r < 2z, but singular when r > 2z, reflecting the singu- 
larity at T — 2z inherent in the correlation functions 
(Ei(x,t')Ei(x,t")) R given in Eqs.(4) and (5). 

For the present purpose, it suffices to show only the 
result of {Av z 2 } for r > 2z [I], 

(Av z 2 ) (6) 

e 2 (t. {t + 2z\ 2 8(1-^) ^ .1 
= 32^p ln {7^2z) + ^^ + 0{p) f 

- - f \ „ (l + -) (for r > 2z) , (7) 
47r 2 m 2 2; 2 V PJ 

where p (> 0) is a dimension- free asymptotic parameter 
for handling the singular integral properly (see the next 
section for details). Accordingly the above expression 



should be understood as an asymptotic expression as p ~ 
0. 

This result is derived by a formula for an asymptotic 
principal value, the rigorous definition of which shall be 
given in the next section; 

P(P) l ^ (X 2 - CT 2 ) 2 

8cr J \1 — a J 2o l p 

for < a < 1. One can derive Eq.(8) with the help 
of Theorem 1; a more direct derivation is also found in 
Appendix C of Ref. [1] . 

Leaving the rigorous treatment of singular integrals 
for the next section, we here focus on the physical rea- 
son why the singularity of correlation functions occurs 
at \T\ = 2z. Due to the mirror-reflections of signals 
with the light velocity, the values of the electric field 
at the two world-points (t',x,y,z) and (t",x,y,z) are 
expected to be strongly correlated when \t' — t"\ = 2z. 
When the measuring time r is short enough (shorter than 
the travel-time of the signal 2z), then it always follows 
\t' — t"\ < 2z, so that these correlations are not captured 
by the probe. When the measuring time is long enough 
(r > 2z), however, these strong correlations accumulate 
in the velocity fluctuations of the particle at z. There- 
fore it is expected that the resulting singular term of the 
form A/p (A > 0) contains information on the reflecting 
boundary. 

On the other hand, typical regularization procedures [5] 
correspond to discarding such a singular term (e.g. the 
1/p term in Eq.(7)) in effect. It should be clarified when 
this type of regularization is valid and when not. We 
shall discuss on this point in more detail in Sec. IV. 

It turns out that the model given here is too simplified 
and should be modified taking into account the switch- 
ing effect [I] and the smearing effect due to the quantum 
spread of the probe-particle [2]. However, it suffices for 
the present purpose of giving some example of singular 
integrals. 

III. BASIC FORMULAS FOR HANDLING 
SINGULAR INTEGRALS 

In view of the example in the previous section, it is 
clear that we sometimes need to estimate a singular in- 
tegral whose integrand possesses a higher-order pole. In 
order to investigate various regularization methods later, 
we first need some concrete quantity corresponding to a 
singular integral for which all the information is retained 
and nothing is discarded. Then, the following asymptotic 
definition of a singular integral may be relevant. 

Definition 1 

Let f(x) be an arbitrary real function defined around 
an interval [a, b], diffcrcntiable at x = c (a < c < b) 
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sufficiently many times. For a positive integer n, then, 
let us introduce an asymptotic principal value of order p 
defined by 



[ b fix) 
P(c,p)U,n) := p( p ) / - tt~ ax 



c-p rb 

Jc+p j 



, f{x \ dx , (9) 
{x-c) n K ' 



where p is a sufficiently small positive parameter. <0 

The asymptotic principal value is a generalization of 
the standard Cauchy principal value, corresponding to 
lim p j,o P(c,p)(/; 1); m two ways. First, the order of singu- 
larity n can be greater than 1. Second, only the asymp- 
totic behavior as p ~ is concerned and the convergence 
for the limit p — > is not necessarily required. In other 
words, we focus on how the integral behaves near p ~ 
rather than the p — > limit itself. In this sense, all the 
information is retained and no infinities are discarded in 
defining the asymptotic principal value. 

Let us now introduce another asymptotic quantity: 

Definition 2 

With the same premises as in Definition 1, we define 

P(p) 





b 


\ m ' 


c-p 

+ 


\ m ] 


_{x - c) n _ 


a 


(x - c) n _ 


a 


(x - c) n _ 



c+p 



It is easily shown that 

Pip) 



\ m - 


b 


\ m i 


_{x - c) n _ 


a 


(x - c) n _ 



{/}& p) ( 10 ) 



with 



{f}{n) } ■= f(c + P )-(-) n f(c-p) . (11) 



By a Taylor-expansion in p. it is obvious that {f}\^^ 
is O(p) for an even n (when /'(c) ^ 0), so that the term 
^{/}(n) P) in Ec l-( 10 ) is Oi 1 /^ 1 ) and singular in the 
p — y limit. Similarly, {/}[^ P ' ) is 0(1) for an odd n 

(when /(c) ^ 0), so that j^{f} ( ^ ] is 0{l/p n ) and sin- 
gular as p — > 0. 

The following definition is just for making formulas 
below concise. 

Definition 3 



C(c,p)(fi 



P{ P ) 



(x - c) ' 



(12) 



With these preparations, let us start with the following 
lemma: 



Lemma 1 



For any function f(x) diffcrcntiablc at x = c and for 
any integer n (> 2), it follows that 

P(c,p)(f,n) = ——^P(c, P ){f\n- 1) - C(c,p)(/)«- 1) • 

(13) 



Proof : 



Noting that 



(x — c) n n — 1 \ (i — c) n 1 



we have 



P(c, P )(f,n) := -^JP( P )J a ( (a._c)n-i ) dx 



n - 1 





[ f(x) 1 


jp(p) 


(x - c) n ~ 1 _ 



-p(p) 



b f i 



a (X-C)"" 1 



dx 



(14) 



where the partial-integral has been performed to get the 
last line. Then the equality follows. □ 

We now prove a formula which relates a multi-pole 
integral with a simple-pole integral: 



Theorem 1 



For any function f{x) differentiable sufficiently many 
times at x = c and for any integer n (> 2), it follows that 



P(cp)(f,n) 
1 



^P(c„)(/^,1) 

-Et^w/^U-*) . as) 

fe=l 1 >' 



Proof : 



(1°) For n = 2, the claimed equality reduces to 
P(c,p) {fi 2) = p( c ,p) (/', 1) - C(o,p) (/> 1) , 

which clearly holds due to Lemma 1. 

(2°) Let us assume that the equality holds for a func- 



tion F(x) and for n = to (m > 2), i.e., 

P(c, P ) {F, m) 

TO — 1 / I SI 

-EE4^)( F(M ' m -^ ■ (16) 

Now applying Lemma 1 for n = to + 1, we have 
P{c, P ){f,rn + 1) = — p( c , p ){f,m) - C(c,p) (/, H 
to ^ (m — lj! v ' 

" £ f^lf C(c) (/ W , m - fc) | - C(c,) (/,"»). 

where the assumed equation (16) for ^(x) = f'(x) has 
been used to get the second equality. Rearranging the 
summation, the last equality reduces to 

P(c,p)(f,m+ 1) 

= 4m(/ w ,i) 



times at x = c and for any integer n (> 2), it follows that 

P(c, P )(f,n) = p( c ,p)(/,n) 

n — 1 

(n-fc-l)l 1 ft _ ni ( c , rt , 



fe=i 



Proof : 



It is straightforward to show this formula due to The- 
orem 1 along with Eqs.(12) and (10). □ 



Lemma 2 



For any function f(x) differcntiable at x = c and for a 
positive integer n, it follows that 



ml 



_^ ((m+l)-fc)! c ^ (/(fc _ 1))(m + 1) _ fc) a cP(c , P )(/,n)=9, 



, P ) (/, n + 1) = - d c P(c . p) (/, n) + . (18) 



Proof : 

We compute d c p( c , P )(f, n) directly as 

c— p />b I 



k=l 



Thus the claimed equation (15) holds for n = to + 1. 

(3°) By the mathematical induction, Eq.(15) holds for 
any integer n (n > 2). □ 

Based on Theorem 1, it is natural to introduce a quan- 
tity p( c ,p) {f, n), which is a simple-pole part plus a regular 
part of p(c,p)(f,n), putting aside singular contributions 
from higher-order poles: 

Definition 4 

We define 



P{e, P ){f,n) := r n _ 1 M p(c, P )(f {n V) 



a */ c+p 



/(*) 



(a: - c) n 



da; 



1 



= np (CjP )(/,n + l)--{/}£ c ^ , 

where the second term in the last line comes from the 
c-derivative applied to the upper- and the lower-limit of 
the integral region. Thus the claimed equation follows. 
□ 



Theorem 3 



For any function f{x) differcntiable sufficiently many 
times at x = c and for an integer n (> 2), it follows that 



1 



P(cp) (/, n) = ^— — dc 11 p (c , p) (/, 1) 



(19) 



{n-k- 1)! 

^ (n-l)\ 

fc=i v ' 



f k - l) {x) 



(x — c) r 



The "mild part" p( c . P )(/, n) of p( c ,p)(/,") shall be im- 
portant in the discussion of regularization methods in 
Sec. IV below. 

We now have a formula which enables us to separate 
singular contributions from a multi-pole integral: 



Theorem 2 



For any function f{x) differcntiable sufficiently many 



Proof : 

Due to Theorem 2, the claimed equation (19) is equiv- 
alent to 

P(c, P ){f,n) = - j-rj d c n ^ p( c ,p)(f, 1) 

^ (n-fc-1)! 1 r f (k-ih M (2Q) 
k=l y >• r 

Thus it suffices to show Eq.(20). 



5 



(1°) Let us consider the case n = 2, where the R.H.S. 
(right-hand side) of Eq.(20) becomes 



d_ 

dc 



dx 



c+p 



+-{f(c-p)+f(c + p)} 

p 



Let us consider a typical singular integral 

/(*) 



I 



a (X-C)» 



dx 



(22) 



for any function f(x) differentiable sufhciently many 
times at x = c (a < c < b) and for a positive integer 
n. 



In this expression, the c-derivative applied to the uppcr- 
and the lower-limit of the integral region yields a term 
which exactly cancels the second term. As a result, the 
above expression reduces to pr CjP )(f, 2), i.e., the L.H.S. 
(left-hand side) of Eq.(20). Thus Eq.(20) holds for n = 2. 

(2°) Let us now assume that Eq.(20) holds for n = m 
(m > 2), i.e., 



P(c, P ) (/, m) = r— [)! dc m P(c,p) (/, 1) 

"f^ (rn-k-l)\ 1 (fc _ x) M 

+ 2^ fm-lV ™-fc u J(m-fc) ^ Z1 J 
fe=i v ; ' r 

Due to Lemma 2, then, it becomes 

P(c,p)(/,™ + l) = — d c p {c , p ){f,m) + — ^{/}^' p) 



m 

= -d, f 

m 



mp" 



'y- 1 (m-fc-1)! 1 r .( fc -i),(c,p) 1 
( m -l)! pm-kU >(m-k)( 



mp" 



■{/} 



where Eq.(21) has been used to get the last line. Noting 
that the relation 

dc{f}£ P) = {/'}ln' P) , 

which obviously holds from Eq.(ll), it reduces to 
P(c, P )(f,m+l) = — d c m p( CjP )(/,l) 



^ (m-fc)! 1 ff ffc-ih( 



fe=l 



c,p) 

m-fc+1) 



Thus Eq.(20) holds for n = 77;, + 1. 

(3°) By the mathematical induction, Eq.(20) holds for 
any integer n (n > 2). Thus the claimed formula Eq.(19) 
has been shown. □ 



IV. TYPICAL REGULARIZATION METHODS 
AND THEIR MUTUAL RELATIONS 

Based on the results in the previous section, let us now 
come back to the problem of regularization methods for 
singular integrals. 



A. Regularization method with partial integrals 

The first method of regularization we consider is a 
method with partial integrals which is sometimes made 
used of. We insert an identity 



1 



1 



1 



(x — c) n n — 1 \ (x — c) n 1 
into Eq.(22) and formally perform a partial integral: 

1 t b ( 1 



n- 1 



(x-c) 
n — 1 ] J a (x — c)" _1 



n-l 



1 



dx 



/(x) dx 

m 

(x - c)™" 1 



(23) 



In this way, the order of singularity is reduced by one. 
Repeating the similar procedure, the integral / is reduced 
to the n = 1 case for which the prescription of the Cauchy 
principal value may be applied. 

Due to Theorem 2, however, it is obvious that singular 
terms should exist and should have been discarded by 
hand in the above procedure. Indeed, compared with 
Eq.(23) with the rigorous expression Eq.(14), it is obvious 
that the singularities which should reside in the second 
term on the R.H.S. of Eq.(23) are simply discarded by 
hand. Thus, in view of Theorem 2, the above method is 
equivalent to the simple replacement of I as 



I ^ 



P(c, P )(f,n) 



(24) 



There is still room, however, to regard the method 
of partial integrals as a shorthand prescription of what 
we here call the method of infinitesimal imaginary part, 
which is much more of theoretical grounds [5]. We shall 
consider this method in the next subsection. 



B. Regularization method with infinitesimal 
imaginary part 

The method of infinitesimal imaginary part is based 
on well-known Dirac's formula [6] for an integral kernel, 



1 



(x — c) ± ip 



P(c,p) 



1 



x — c 



=F iir5(x — c) 



(25) 
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which is most easily shown by estimating an integral 
la (x-c)±ip ^ means of an appropriate contour-integral 
for a suitable function f(x). 

By differentiating the both-sides of Eq.(25) n — 1 times 
with respect to c, and by applying Theorem 3, we get 



— > 2 in the limit p — > 0. This is because all the infor- 
mation is retained in the prescription of the asymptotic 
principal value. 

On the other hand, the above mentioned scheme makes 
a replacement 



{{x-c)±ip) n 



(n-1)! 



S^-^ix-c) (26) 



in the sense of an integral kernel. Recalling Definition 4, 
however, we see that the R.H.S. is reduced to the n = 1 
case, for which the prescription of the Cauchy principal 
value may be applied. 

It is notable that just the introduction of some in- 
finitesimal imaginary part results in a tamable quantity 
such as p(p. c ) ( ■ , n) at the cost of the imaginary con- 
tribution of the second term on the R.H.S. of Eq.(26). 
Thus along with some causality arguments [5] , it is often 
argued that the singular integral / should be interpreted 



as the real part of f 



{{x-c)±ip) 



:dx, i.e., 



/(*) 



((x - c) ± ip) 



rdx 



(27) 



As far as one is evaluating real quantities, one may fur- 
ther argue that the second term on the R.H.S. of Eq.(26) 
shall not contribute. If so, the procedure is in effect 
equivalent to the replacement Eq.(24). In this sense, the 
method of partial integrals discussed in the previous sub- 
section may be justified provided that it is regarded as 
a shorthand prescription of the method of infinitesimal 
imaginary part. 

Another way of looking at this method is to pay atten- 
tion to the L.H.S. (rather than the R.H.S.) of Eq.(26). 
As far as computations of real quantities are concerned, 
then, this method is equivalent to the replacement , J 1 ,„ 

[X C) 

with ( ^ x ^i +p -i ^ along with taking the limit p — >• af- 
ter evaluating the integral: 



1 1— y lim 



fix) 



■r 



dx 



(28) 



There is some subtle points in this method. One of 
them is to discard the imaginary part of the R.H.S. 
of Eq.(26) on the grounds that one is evaluating real 
quantities. Considering that the regularization has been 
achieved at the cost of introducing the imaginary part 
though tiny, the imaginary part should carry important 
information and some concern naturally arises whether 
one can discard it so freely. 

Indeed, a simple example can be presented for which 
this kind of procedure fails. Let us consider an integral 
1\ = J_j dx which is purposefully regarded as 



h 



x ■ —dx 

l x 



It is obvious that I 1 = 2. The analysis by the asymptotic 
principal value (see the next subsection) also results in 



1 x 

— H-> 



so that 



x x 2 + p 2 



t r(p) — 
i\ — > i l .— 



= I(ln(, 2 + p 2 ))' , 

l i 
x ■ - (\n(x 2 + p 2 ))' dx 



= ln(l + p 2 ) - f ln(x 2 + p 2 )dx , (29) 

J<3 

where a partial integral has been performed to get the 
last line. However, it is clear that I± — > oo as p —> 0, 
contradicting with the obvious result I\ = 2. 

Quite interestingly, no contradiction occurs for I m = 
J j x m ■ ^dx with m > 2 since the second term in Eq.(29) 

becomes — m x m_1 ln(x 2 + p 2 )dx so that no singular- 
ity occurs around x ~ for m > 2. More generally, 
the integral of the form ^^-dx, if treated by the 
above prescription, gives rise to the dominant contribu- 
tion — /'(0) Jq 1 ln(2; 2 + p 2 )dx which diverges as p — > 0. 

With these caveats in mind, let us now move to a new 
regularization method based on the asymptotic principal 
values. 



C. Regularization method with asymptotic 
principal values 

Let us finally introduce a new regularization method 
based on the asymptotic principal values. 

For the simple example in Sec. II, there has been a def- 
inite physical interpretation of the singularity in the cor- 
relation function. Furthermore, the system considered 
there has been a combination of quantum objects with 
a macroscopic mirror. Therefore it might be also prob- 
able that the deepest cause of the singularity resides in 
the validity issue of the model originating from too much 
extrapolation from the quantum side to the macroscopic 
situation. Indeed there is an investigation showing that 
the quantum fluctuations of the mirror boundary drasti- 
cally decrease the singular behavior near the mirror [7]. 
Therefore it is reasonable to take the origin of the sin- 
gularity more realistically (rather than just mathemati- 
cal phenomenon) , expecting that some physical processes 
suppress the order of singularity. 

Going back to the example of the integral I in Eq.(22), 
then, it is possible to interpret / in the sense of an asymp- 
totic principal value, 



I >-> P(c, P )(f,n) 



(30) 



with the dimension-free parameter p being provided by 
the ratio of some natural cut-off scale with the system- 
size in question. (For instance, the ratio of the plasma 
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wave-length of the mirror with 2z for the example in 
Sec. II). The advantage of this regularization scheme is 
that one can explicitly analyze the p-depcndence of the 
integral. For instance, one may study the influence of 
the quantum fluctuations of the mirror by treating p as 
a fluctuation parameter. The result for (Av z 2 ) given in 
Eq.(7) along with Eq.(8) is an example of the computa- 
tion by the method of asymptotic principal values. 

We see that Theorem 2 is the basis for understand- 
ing the relation between the regularization methods dis- 
cussed so far. The difference between the method of 
asymptotic principal value (p( c , p )(/, n)) and the method 
of infinitesimal imaginary part (p( c , P )(/, n)) is given by 
the second term on the R.H.S. of Eq.(17), which is of 

oa/p"- 1 )- 

V. SUMMARY AND DISCUSSIONS 

In this paper, we have focussed on singular integrals 
with a higher-order pole which frequently emerge in com- 
puting quantities based on two-point correlation func- 
tions of a vacuum. 

To deal with this type of singular integrals, we have 
introduced the concept of asymptotic principal values. 
The asymptotic principal value of order p, which is a 
generalization of the Cauchy principal value, is defined 
by introducing a cut-off parameter p, focussing solely on 
the asymptotic behavior of the integral as p ~ 0. In this 
sense, it is a rigorous object retaining all the information 
on the singular integral. 

We have then proved several theorems on asymptotic 
principal values which are expected to serve as bases for 
studying regularization methods for singular integrals. 

To see how asymptotic principal values can be made 
use of, we have selected three typical regularization meth- 
ods and have analyzed their mutual relations with the 
help of theorems we have prepared. It has turned out 
that the concept of asymptotic principal values and re- 
lated theorems are quite useful in this kind of analy- 



sis. Indeed, in terms of asymptotic principal values, it 
has been possible to describe without ambiguity what 
is discarded and what is retained in each regularization 
method. 

No universal regularization method is available so far 
and we need to carefully select or invent a suitable 
method depending on the problem in question. For in- 
stance, we recall the example in Sec. II where velocity 
dispersion of the probe, (Av z 2 ), is sensitive to the reg- 
ularization method. In particular, the result expected 
by the method of infinitesimal imaginary part (Sec. IV B) 
(and the method of partial integrals (Sec. IV A)) is 

(Av z 2 ) ~ f 2 - (forr»2z). 

On the other hand, the result expected by the method of 
asymptotic principal values (Sec. IV C) is 

(AO~ f 2 J l + -) (forr»2z) , 

by choosing p in the order of the ratio of plasma wave- 
length and the typical size 2z. Strictly speaking, the 
model in Sec. II is a too simplified one and should be 
modified taking into account the quantum spread of the 
probe-particle itself. Then the behavior of (Av z 2 ) at late 
time is corrected to a more reasonable one (Av z 2 ) ~ 1/r 2 
rather than ~ 1/z 2 [2]. 

In any case it is significant to compare the results de- 
rived by different regularization methods in more detail 
for approaching to a more satisfactory mathematical the- 
ory of regularization procedures. 
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